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Abstract 

Ferrofluids containing nonmagnetic particles are called inverse ferrofluids. Based on the Ewald-Kornfeld 
formulation and the Maxwell-Garnett theory, we theoretically investigate the magnetophoretic force exert- 
ing on the nonmagnetic particles in inverse ferrofluids due to the presence of a nonuniform magnetic field, 
by taking into account the structural transition and long-range interaction. We numerically demonstrate 
that the force can be adjusted by choosing appropriate lattices, volume fractions, geometric shapes, and 
conductivities of the nonmagnetic particles, as well as frequencies of external magnetic fields. 
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I. INTRODUCTION 



Magnetizable particles experience a force in a nonuniform magnetic field, which comes from 
the interaction between the induced magnetic moment inside the particles and the external mag- 
netic field. This phenomenon is called magnetophoresis/'^ and has been exploited in a variety of 
industrial and commercial process for separation and beneficiation of solids suspended in liquids.^ 
Magnetophoresis may also open up a way to characterize and separate cells for biochemical anal- 
ysis based on intrinsic and extrinsic magnetic properties of biological macromolecules.^ 

Ferrofluids, also known as magnetic fluids, are a colloidal suspension of single-domain ferro- 
magnetic nanoparticles (e.g., made of Fe304), with typical dimension around lOnm, dispersed 
in a carrier liquid like water or kerosene.^'^"^ Ferrofluids containing nonmagnetic particles (or 
magnetic holes) are called inverse ferrofluids. So far, inverse ferrofluids have drawn much atten- 
tion for their potential industrial and biomedical applications. ^^^^ The size of the nonmagnetic 
particles (such as polystyrene particles) is about 1~100 /xm. Since the nonmagnetic particles 
are much larger than the ferromagnetic nanoparticles in host ferrofluids, the host can theoreti- 
cally be treated as a one-component continuum in which the much larger nonmagnetic particles 
are embedded. ^^'^^ It is known that the ground state of electrorheological fluids and magnetorhe- 
ological fluids is body-centered tetragonal (bet) lattices. ^^'^^ Due to mathematical analogy, the 
ground state of inverse ferrofluids might also be body-centered tetragonal (bet) lattices. Similar to 
electrorheological fluids, under appropriate conditions, such as by adjusting the application of 
external magnetic fields, structural transition may also appear in inverse ferrofluids, e.g., from the 
bet to the body-centered cubic (bcc), and then to the face-centered cubic (fee) lattices. 

When an external magnetic field is applied, the nonmagnetic particles in inverse ferrofluids can 
be seen to possess an effective magnetic moment, which, however, is anti-parallel to the magneti- 
zation of the host ferrofluid. In this sense, the nonmagnetic particles also experience an equivalent 
magnetophoretic force in the presence of a nonuniform magnetic field naturally. In this work, we 
shall apply a nonuniform magnetic field to investigate the magnetophoresis of the nonmagnetic 
particles in inverse ferrofluids with different structures of specific lattices. 

One of us^^ has investigated the dielectrophoretic force acting on a microparticle in electrorheo- 
logical solids under the application of a nonuniform external electric field by using a self -consistent 
method^^ and the spectral representation approach. In our present work, we shall use the Ewald- 
Kornfeld formulation and the Maxwell-Garnett theory for uniaxially anisotropic suspensions to 
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calculate the effective permeability of the inverse ferrofluid, taking into account the long range 
interaction arising from specific lattices. Based on this, the magnetophoretic force acting on the 
nonmagnetic particles can be calculated accordingly. We shall focus on the role of types of lat- 
tices, volume fractions, geometric shapes, and conductivities of the nonmagnetic particles, as well 
as frequencies of external fields. 

This paper is organized as follows. In Section UIl we first briefly review the theory of mag- 
netophoresis (Section FlIAl) . By using the Ewald-Komfeld formulation and the Maxwell- Gamett 
theory for uniaxially anisotropic suspensions, we calculate the effective permeability of the in- 
verse ferrofluids (Section Til B|) . Then we investigate magnetophoresis with eddy current induction 
(Section HTCl ). In SectionHIIl we give the numerical results under different conditions. This paper 
ends with a discussion and conclusions in Section HVl 

n. FORMALISM 

A. Theory of magnetophoresis 

To establish our theory, we first briefly review the theory of magnetophoresis (for details of 
magnetophoresis, we refer the reader to ref 1). Let us start by considering a homogeneous sphere 
with radius R and magnetization Mi that is suspended in a magnetically linear fluid of permeabil- 
ity H2 and subjected to a uniform magnetic field Hq. The relation between magnetic flux density 
B, magnetization M, and magnetic intensity H is given by 



where /io is the permeability of free space. Eq \T\ is completely general, the magnetization M 
includes permanent magnetization plus any linear or nonlinear function of H. Mi is assumed to 
be parallel to Hq, Mi ||Ho. We also assume that no electric current flows anywhere. Thus, there is 



everywhere. Then, this magnetostatics problem may be solved using a scalar potential ip, defined 
by H = —Vip. The general solutions for 'ipi and ip2 outside and inside the sphere are given by 
respectively, 



B = /io(H + M), 



(1) 



V X H = 



(2) 



Mr, 9) 




r > R. 



(3) 



br cos 9, r < R. 



(4) 
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The coefficients a and b can be determined by considering the appropriate boundary conditions, 

„_ f^-f'.R^Ho + —^R'M,, (5) 



Ho + 2fi2 /iO + 2/i2 

b=^^Ho-—^M,. (6) 

/iO + 2/i2 /iO + 2/i2 

It is worth noting that b denotes the magnitude of the uniform magnetic field Hi inside the sphere. 
We can determine the effective magnetic dipole moment niefr by comparing the dipole term in eq|3] 
with 



rriefi cos 9 
Anr 

Then, the effective magnetic moment is given by 



(8) 



jiQ + 2/i2 A^o + 2/i2 
For a magnetically linear particle of radius R and permeability fii suspended in a host medium 
with permeability yU2, we have 

Ml = xiHi, (9) 

where xi = /^i/z^o — 1 is the susceptibility of the particle. In this case, the effective magnetic 
dipole moment vector niefr takes a form very similar to that for the effective dipole moment of 
dielectric spheres, 

meff = 47ri?3ir(/ii,/i2)Ho, (10) 
where the Clausius-Mossotti factor K^hi, 112) is given by, 

K{HU^2) = -^^^^. (11) 

Ail + 2/i2 

The magnetophoretic force exerted by a nonuniform magnetic field on a magnetic dipole can then 
be written as 

i^dipoie = At2mcfr ■ VHq. (12) 

Combining this expression with that for the effective magnetic dipole moment eq [lOl the mag- 
netophoretic force exerted on a magnetizable spherical particle in a nonuniform magnetic field is 
written as^ 

Fmap = 27r/i2i?'i^(/ii, ^^2)VHl. (13) 

Eq[T3]is similar to the expression of the dielectrophoretic force on dielectric particles. ^'^^ It is 
apparent to see that the magnetophoretic force are proportional to the particle volume and the 



permeability of the suspension medium. And the force is directed along the gradient of the external 
magnetic field. There are two kinds of magnetophoresis depending upon the relative magnitudes of 
/ii and /i2. If K > (namely, fii > 112), particles are attracted to magnetic field intensity maxima 
and repelled from minima. This phenomenon is called positive magnetophoresis. On the contrary, 
if K < (i.e.. Ill < 112), inverse behavior appears, which is called negative magnetophoresis. 
Generally, there may be two magnetophoresis systems. One is that ferromagnetic particles are 
suspended in a medium that can be either nonmagnetic or less (linear) magnetizable, and the other 
one is that nonmagnetic particles are in a magnetizable liquid like ferrofluids. As for the second 
case of our interest in this work, in case of a moderate magnetic field (satisfying the assumption of 
a magnetically linear host), we can use eq[T3]by setting /ii = /xq to calculate the magnetophoretic 
force on the nonmagnetic particles submerged in a host ferrofluid. 

B. Effective permeability 

We resort to the Ewald-Komfeld formulation to calculate the effective permeability by includ- 
ing the structural transition and corresponding long-range interaction explicitly. A bet lattice can 
be regarded as a tetragonal lattice, plus a basis of two particles each of which is fixed with an 
induced point magnetic dipole at its center. One of the two particles is located at a comer and the 
other one at the body center of the tetragonal unit cell. The tetragonal lattice has a lattice constant 
ci= along the z axis and lattice constants ai = hi = ^q^^ along x and y axes, respectively. 
In this paper, we identify a different structure transformation from bet to bcc, and to fee lattices 
by changing the uniaxial lattice constant ci under the hard-sphere constraint. The volume of the 
tetragonal unit cell Vc remains unchanged (i. e., Vc = ^^) as q varies. In this way, the degree of 
anisotropy of the tetragonal lattices is measured by how q is deviated from unity and the uniaxial 
anisotropic axis is along the z axis. In particular, q = 0.87358, 1, and 2^ represent the bet, bcc, 
and fee lattices, respectively. 

When an external magnetic field Hq is directed along the x axis, the induced magnetic moment 
iTiefr is perpendicular to the uniaxial anisotropic axis (z axis). At the lattice point R = 0, the local 
field can be determined by using the Ewald-Komfeld formulation,^°~^^ 

= meff^E[-7i(^.)+^,V72(^i)] 

j=l R^O 
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^^i:n(G)iexpfz^Ui^. (14) 



where 71 and 72 are two coefficients, given by 

, , erfcirir) 2ri / 2 2\ 



3erfc{r]r) ^rf 6r] ^ ^ 

where erfc(?7r) is the complementary error function, and r] is an adjustable parameter making the 
summations in real and reciprocal lattices converge rapidly. In eq [l4l R and G denote the lattice 
vector and reciprocal lattice vector, respectively, 

where I, m, n, m, v, and w are integers. In addition, Xj and Rj in eq[l4]are given by 

j - 1 

= |R —{ai± + aiy + ciz)\, 

andn(G) = l+exp[i{u+v+w) /n] is the structure factor. The local magnetic field can be obtained 
by summing over all integer indices, {u,v,w) 7^ (0,0,0) for the summation in the reciprocal 
lattice, and (j,l,m,n) 7^ (1,0,0,0) for the summation in the real lattice. In our calculation, 
because of the exponential factors, we may impose an upper limit to the indices, i.e., all indices 
ranging from —L to L, where L is a positive integer. We consider an infinite lattice. So, for finite 
lattices, one must be careful about the effects of different boundary conditions. 
Now, let us define a local field factor a as 

3 Hl 3V,Hl .... 

where J2y^ nieff = 2meff means the sum of the two vector magnetic moments in the unit cell. 
It is worth remaking that a is a function of single variable q, and that it stands for ar (local 
field factor in transverse field cases) and aL (local field factor in longitudinal field cases) which 
satisfy a famous sum rule in condensed matter physics for ar and aL, 2aT + ai = 3.^^'^^ Here 
the longitudinal (transverse) field case corresponds to the fact that the external magnetic field 
is parallel (perpendicular) to the uniaxial anisotropic z axis. Next, for calculating the effective 
permeability, we use the Maxwell-Gamett theory for uniaxially anisotropic suspensions, ^^'^^ 

/^e — ^ /ii — /i2 ^^^^ 



(^L,Tl^e + (3 - 0;L,T)/i2 /^l + 2/i: 
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where yUg stands for the effective permeability of the inverse ferrofluid, /xi the permeability of 
the nonmagnetic particles, namely fii = yUo, Ai2 the permeability of the host ferrofluid, and / the 
volume fraction of the nonmagnetic particles in the inverse ferrofluids. 

It should be noted that eq [16] holds for the case of nonmagnetic particles with a spherical shape 
only. For nonspherical cases (e.g., prolate or oblate spheroids), we may use demagnetizing factor 
g to describe the shape of different kinds of spheroids. For a prolate spheroid with three principal 
axes, a, b, and c, there is a = b < c, and the demagnetizing factor along the major axis c is^^'^^ 



1 



gL 



k 



:\n{k + vW 



(17) 



where k = c/a > 1 and g^ < 1/3. For an oblate spheroid with a < b = c, the demagnetizing 
factor g^ along the major axis a is^^'^^ 



9l 



1 - 



■ arcsm ■ 



(18) 



.(A;2-l)t k k^-l_ 

where k = c/a > 1 and gL > 1/3. For both cases, there is a geometrical sum rule for gL and 
gr, 2gT + (^l = 1, where gx denotes the demagnetizing factor along the minor axis. Note that the 
spherical shape of nonmagnetic particles corresponds to (7/, = (7^ = 1/3. So far, to include the 
nonspherical shape effects, erq[T6]for calculating the effective permeability of inverse ferrofluids 
containing spherical nonmagnetic particles might be modified as 



Apparently, the substitution of gL = gr = ^/'^ into eq[T9] reduces to eq[T6l 



/- 



(19) 



C. Magnetophoresis with eddy current induction 

When an ac (alternating current) magnetic field is applied to an electrically conductive magne- 
tizable particle, eddy currents are induced, which tend to oppose penetration of the field into the 
particle. At sufficiently high frequencies, these eddy currents strongly influence the effective mag- 
netic moment. Within the particle, V x H 7^ 0. Similar derivation of this section can be found in 
ref 1 for a magnetic particle in a linear medium. Nevertheless, here we discuss a different system, 
thus yielding the following formulae with some difference. 

Let us consider a nonmagnetic conductive particle with an effective magnetic dipole of magni- 
tude nies [note that in this case there is nieff || (— Hq) as mentioned in Section[ll, which is embedded 
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in the system containing both a magnetically linear host ferrofluid and many other nonmagnetic 
conductive particles. In this case, the nonmagnetic particles are also assumed to situate in specific 
lattices and exist in the shape of spheres. The effective permeability of the whole system /ig can 
be calculated according to eq[T6l The magnetic vector potential ^ at position r may be written as 

* = /\ , (20) 

where V x ^' = B, B is the magnetic flux density. If we align the moment anti-parallel with z 
axis, that is mefr = —niesi- Then eq|20]can be rewritten as 

LLemeS sin 9 7 ^^^^ 
* = - w (21) 

where (f) is the azimuthal unit vector in spherical coordinates. We consider a conductive nonmag- 
netic spherical particle of permeability /Ui(= /ig) and conductivity cti, which is subjected to a 
uniform, linearly polarized magnetic field H(t) = Hq cos cotz. The effective moment of the sphere 
can be written as 

mcff = -A7iR^L{uj)Hoz, (22) 

where L{uj) = R^D /2 is a frequency-dependent magnetization coefficient, and D depends on /Xe, 
111, ai, and uj, 

^ ^ (2/il + fie)^I-l/2 - [/ie(l + ^'^) + 2/Xl]/+i/2 ^^3) 
(/il - He)l^I-l/2 + K(l + ^'^) - l^l]I+l/2 



with I±i/2 = I±(u) and u = j^iaiOjR. The quantities /±i/2 are half-integer order, modified 
Bessel functions of the first kind. For a nonuniform ac magnetic field, according to the effective 
dipole method (see Sections III AI and IIIB|) . the time-average magnetophoretic force acting on a 
nonmagnetic conductive particle is 

(FMAp(t)) = yReKff ■ VH]. (24) 

Next, we can obtain the expression for the magnetophoretic force in terms of the real part of L, 

{FuApit)) = -27rfi,R'Re[L]VHl (25) 

A simplified expression for Re[L] may be obtained by invoking certain Bessel function identities, 

^ 3fiRifiR-l)B' + Rl{smh2Rm-sm2Rm)A 1 

^'^^^ = ~2 im-iyB^ + RlA- 2' ^^^^ 
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where fiR = /ii//ie is the relative permeability, Rm = -Ryco'/^iCi/S a dimensionless modulus 
called magnetic Reynolds number, and A and B are respectively given by 

A = 2i?m(cosli 2i?m — cos2i?m) + (yU/j — l)(sinh2_Rm — sin2_Rm), 
B = cosh 2Rm — cos 2Rm — -Rm(sinh 2Rm + sin 2Rm)- 

III. RESULTS 

If the shape of the nonmagnetic particle is spherical, the magnetophoretic force on the nonmag- 
netic particle Fmap can be calculated by 

Fmap = 27ifi,R'Kifii,fi,)VHl (27) 
Kinuf^e) = ^^if^. (28) 



For obtaining eq |27l in order to include the lattice effect, fi2 in eq \\3\ has been replaced by fie 
that can be calculated according to eq [HI During the numerical calculation, for convenience, 
F stands for — Fmap/ (271 R'^VH^). In our calculation, we choose the parameters /ii = /iq and 
fi2 = 2.9/io- In this case, the parameter /^(yUi, /ig) < corresponds to the phenomenon of negative 
magnetophoresis . 

Figure \\\ displays the dependency of F on g when the external nonuniform magnetic field 
is parallel to the uniaxial anisotropic axis (longitudinal field cases) for three different volume 
fractions / = 0.15, / = 0.25, and / = 0.35. While Figure[2]shows the dependency of F on g when 
the external nonuniform magnetic field is perpendicular to the uniaxial anisotropic axis (transverse 
field cases) for the same three different volume fractions. Because of the existence of anisotropy. 
Figure [U displays that F decreases as q increases. Inverse behavior appears in Figure [2l It is 
apparent that the lattice effect plays an important role in magnetophoresis of nonmagnetic particles 
in inverse ferrofluids. In the mean time, increasing volume fraction / of nonmagnetic particles 
leads to decreasing normalized magnetophoretic forces for both longitudinal and transverse field 
cases. As the volume fraction of the nonmagnetic particles increases, the lattice effect becomes 
stronger. By using the Maxwell- Gamett theory for uniaxially anisotropic suspensions eq \T6[ a 
smaller effective permeability of the inverse ferrofluid can be obtained. Consequently, the absolute 
value of the normalized magnetophoretic force on the nonmagnetic particles becomes smaller. 

The inclusion of the shape effects makes the magnetophoretic force expression more compli- 
cated for a nonspherical particle than for a spherical. For prolate spheroidal particles (a = b < c). 



which are aligned with its major axis being parallel [corresponding to the case of a/, in eq \T9\ 
and (7l in eqs \T9\ and [29l or perpendicular [reflected by ax in eq \T9\ and gx in eqs \T9\ and [2911 
to the direction of the external nonuniform magnetic field, the magnetophoretic force expression 
resembles that for a sphere, 

VHl (29) 

where fie should be calculated according to eq [19] by including the lattice effects. From this 
equation, we find that the force remains proportional to the scalar product of particle volume and 
excess permeability. Ellipsoidal particles exhibit the same positive and negative magnetophoresis 
phenomenon, depending on whether fii is greater than or less than /ig- This expression eq[29|is 
also valid for oblate spheroids (a < 6 = c) by replacing a^c with ac^. In fact, a particle with its 
short axis (e.g., a axis for the oblate spheroid discussed herein) parallel to the magnetic field is 
in a very unstable alignment. In this work, for completeness of including shape effects, we shall 
also discuss such a case of oblate spheroidal particles with its short axis (i.e., major axis or a axis) 
parallel to the magnetic field. For a spherical particle, there is a = b = c and gL = gr = ^/'^, and 
eq[29]reduces to eq [27] naturally. 

Next, for convenience, we denote Fi as — FMAp/(27ra^cVi/o) for prolate spheroids or 
— -^MAp/(27rac^ViJo) oblate spheroids, which means that we only calculate the normalized 
magnetophoretic force. 

Figure l3] displays the dependency of Fi on q when the external nonuniform magnetic field 
is parallel to the uniaxial anisotropic axis (longitudinal field cases) for different aspect ratios k 
for prolate spheroidal cases. The normalized magnetophoretic force of the nonmagnetic particle 
decrease as g (g > 0.7) or /c increases. Nevertheless, at small q, 0.6 < q < 0.7, the situation is 
complicated, but the difference for different k is very small. On the other hand, for the transverse 
field cases shown in Figure ]4l as q and/or k increases, the normalized magnetophoretic force is 
caused to increase for the full range of q. 

Figure 15] displays the dependency of Fi on q when the external nonuniform magnetic field is 
perpendicular to the uniaxial anisotropic axis (transverse field cases) for different aspect ratios k 
for oblate spheroidal cases. The normalized magnetophoretic force of the nonmagnetic particle 
decreases as q increases, and increases as the aspect ratio k increases. Inverse behavior appears 
for the transverse field cases, which is plotted in Figure[6] In addition, the relationship between Fi 
and the volume fraction of the nonmagnetic particles is also investigated (no figures shown here). 
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-^MAP 



27ra^c 



1 (/^l - f^e) 

31 + 



9l,t 



The normalized magnetophoretic force on the nonmagnetic particle is found to decrease as the 
volume fraction increases, which are consistent with the result of the spherical case. 

In section III CI we consider the magnetophoresis with eddy current induction. The conduc- 
tivity of the nonmagnetic particle (such as polystyrene) can be changed by using the method 
of doping. Again, for convenience, F2 stands for the normalized magnetophoretic force, F2 = 
— (-FMAp(^))/27r_R^V//o. Figure |7] displays the dependency of F2 on different frequencies cu for 
different conductivities a. It is shown that the normalized magnetophoretic force of the nonmag- 
netic particle decrease as uj and/or a increase for the current parameters in use. 

IV. DISCUSSION AND CONCLUSION 

We have theoretically investigated the magnetophoretic force on the nonmagnetic particles sus- 
pended in a host ferrofluid due to nonuniform magnetic fields, by taking into account the effects 
of structural transition and corresponding long-range interaction. We have resorted to the Ewald- 
Kornfeld formulation and the Maxwell-Garnett theory for calculating the effective permeability of 
the inverse ferrofluids. And we have also investigated the role of the volume fraction, geometric 
shape, and conductivity of the nonmagnetic particles, as well as frequencies of external fields. 

In this work, we have plotted the figures by using the values for the volume fraction, / = 
0.15,0.25 and 0.35, which actually corresponds to low concentration systems. Nevertheless, it 
would be realizable as the nonmagnetic particles can have a solid hard core and a relative soft 
coating (e.g., by attaching long polymer chains, etc.) to avoid aggregation. Thus the model of a 
soft particle with an embedded point dipole has been used throughout this work. In this case, the 
many-body (local-field) effects are of importance, while the multipolar effects can be neglected. 
This is exactly part of the emphasis in this work. As the nonmagnetic particles are located closely, 
the multipolar interaction between them is expected to play a role.^'''^^ In this regard, it is of value 
to extend the present work to investigate the effect of multipolar interaction. In addition, it is also 
interesting to see what happens if the particle is superconductor. In the Appendix, we develop 
the theory for describing the magnetophoretic force acting on type I superconducting particles. 
In all our figures, only the normalized magnetophoretic forces are displayed because it is very 
complicated to calculate the field gradient Vi^o the lattices of our discussion accurately. In 
so doing, the contribution of the gradient has not been reflected in the figures. Fortunately, the 
difference of the gradient for the different lattices can be very small due to the closely packing of 
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the nonmagnetic particles in the lattices. 

The present theory is valid for magnetically linear media, e.g., ferrofluids in a moderate mag- 
netic field. Throughout the paper, for convenience we have assumed the host ferrofluid to behave 
as an isotropic magnetic medium (continuum) even though a magnetic field is applied. In fact, an 
external magnetic field can also induce a host ferrofluid to be anisotropic, and hence the degree 
of anisotropy in the host depends on the strength of the external field accordingly. Also, in Sec- 
tion III CI we have assumed that only the suspended nonmagnetic particles are conducting while 
the host is not. Therefore, it is also interesting to investigate the cases of an anisotropic and/or 
conducting host. 

For magnetic materials, nonlinear responses can occur for a high magnetic field. In case of 
moderate field, linear responses are expected to dominate. The latter just corresponds to the case 
discussed in this work. We have also investigated magnetophoresis with eddy current induction. In 
fact, conductive particles can exhibit diamagnetic responses in ac fields because of induced eddy 
currents. Nevertheless, the magnetic susceptibilities of some conductive (diamagnetic) particles 
like silver or copper are generally —10^^ or so. Thus, their relative magnetic permeabilities could 
be 1 + (—10^^) ~ 1, as already used for Fig.|7l 

For treating practical application, the present lattice theoretical model should be replaced by the 
random model, mainly due to the existence of a random distribution of particles. In so doing, some 
effective medium approximations may be used instead, e.g., the Maxwell-Garnett approximation, 
the Bruggeman approximation, and so on. 

To sum up, our results have shown that the magnetophoretic force acting on a nonmagnetic 
particle suspended in a host ferrofluid due to the existence of a nonuniform magnetic field can 
be adjusted significantly by choosing appropriate lattices, volume fractions, geometric shapes and 
conductivities of the nonmagnetic particles, as well as frequencies of external fields. 
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Appendix: Magnetophoresis of type I superconducting particles 

A conventional type 1 superconducting particle exhibiting the Meissner effect may be modelled 
as a perfect diamagnetic object. Such a model is adequate as long as the magnitude of a magnetic 
field does not exceed the critical field strength of the material. If these limits are not exceeded, 
we may estimate the magnetic moment nicfr for a type I superconducting particle by taking the 
K(jj,i, = —0.5 limit of eqs[lO]and[T3l and then obtain the effective dipole moment niefr 

meff = -2rfHo, (30) 

and the magnetophoretic force Fmap 

Fmap = -TifieR'VHl (31) 

Eq [3T] indicates that type I superconducting particles will always seek minima of the external 
applied magnetic field and can be levitated passively in a cusped magnetic field. As a matter of fact, 
eq[3T]is developed from that (see ref 1) for a single type I superconducting particle by taking into 
account the lattices effects or, alternatively, introducing the parameter of effective permeability /ig. 
It is known that the diamagnetic behavior is very strong in a type I superconductor. For such a 
superconducting particle, which is perfectly diamagnetic, its magnetic susceptibility is —1, which 
causes its relative magnetic permeability to be 0. Nevertheless, this does not affect the validity of 
the formulae presented herein. 
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Figure Captions 

Figure 1. Cases of nonmagnetic spherical particles. F as a function of q for different volume 
fractions, for longitudinal field cases. Parameters: /xi = /iq and ^2 = 2.9/io- 
Figure 2. Same as Figured! but for transverse field cases. 

Figure 3. Cases of nonmagnetic prolate spheroidal particles. Fi as a function of q for different 
aspect ratios: = 4, = 6, A; = 8, and k = 10, for longitudinal field cases. Parameters: /ii = fj.Q, 
/i2 = 2.9/^0, and / = 0.15. 

Figure 4. Same as Figure [3l but for transverse field cases. 

Figure 5. Cases of nonmagnetic oblate spheroidal particles. Fi as a function of q for different 
aspect ratios: = 4, = 6, A; = 8, and k = 10, for longitudinal field cases. Parameters: Hi = Hq, 
fi2 = 2.9|Uo, and / = 0.15. 

Figure 6. Same as Figure [51 but for transverse field cases. 

Figure 7. Cases of nonmagnetic, electrically conductive spherical particles. F2 as a function 
of Lu for different conductivities: ai = 10^'^ S/m, 10^^ S/m, and 10^^ S/m. This figure is for bet 
lattices, namely, q = 0.87358 (which yields = 1.09298), and the longitudinal field cases. 
Parameters: /xi = /xq, ^^2 = 2.9/xo, -R = 2.0 x 10"^ m, and / = 0.15. 
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FIG. 1: /Gao, Jian, Zhang, and Huang 
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FIG. 4: /Gao, Jian, Zhang, and Huang 
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FIG. 7: /Gao, Jian, Zhang, and Huang 
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